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Abstract 



We report six-dimensional quantum dynamical calculations of dissociative 
adsorption and associative desorption of the system H2/Pd(100) using an ab 
initio potential energy surface. We focus on rotational effects in the steering 
mechanism, which is responsible for the initial decrease of the sticking prob- 
ability with kinetic energy. In addition, steric effects are briefly discussed. 
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I. INTRODUCTION 



In recent years the dynamics of dissociative adsorption has been the subject of a large 
number of experimental and theoretical investigations (see, e.g., Refs. |1],0J^^]). As far 
as quantum dynamical simulations were concerned, these studies were restricted to low- 
dimensional calculations on model potentials due to computational constraints and the non- 
availability of ab initio potential energy surfaces (PES). By varying potential parameters 
experimental results were tried to be reproduced qualitatively. These studies laid the foun- 
dations of the current understanding of simple surface reactions and of the topological fea- 
tures realistic potentials should have. The main effects of molecular vibration |5|,|£l,[7||^|9lJlO| , 
rotation [|ll|,|l^,|l^,|l^,|l^,|l^] and lateral corrugation [|l7i|l8| on the dissociative adsorption 
probability seemed to be understood to a large extent, but it remained unclear whether 
the qualitative explanations would still be valid in high-dimensional dynamical calculations 
including all crucial degrees of freedom. 

Just recently it has become possible to evaluate the six-dimensional potential energy sur- 
face of hydrogen dissociation on metal surfaces |T9| , |20|j21| , p2[| by density-functional theory. 
This development also enforced new efforts for improving the quantum dynamical algo- 
rithms. Indeed it is now feasible to perform studies of hydrogen dissociation where all six 
degrees of freedom of the hydrogen molecule are treated quantum mechanically . These 
calculations showed that the initial decrease of the sticking probability with kinetic energy 
found experimentally for H2 on Pd(100) |24|] and on many other transition metal surfaces 
p5| , |26|p7| , |28|p9| , |30|| is not due to a precursor mechanism, as was commonly believed, but can 
be explained by dynamical steering. 

In this contribution we will - after briefly recalling the theoretical background and the 
main result of our previous study [^] - focus on the influence of rotations on the adsorption 
dynamics in the system H2/Pd(100). We will describe the dependence of the sticking proba- 
bility on the initial rotational quantum number ji of the impinging hydrogen molecules and 
show how this dependence could be verified experimentally. We end with a brief discussion 
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of steric effects and concluding remarks. 



II. THEORETICAL BACKGROUND 



The potential energy surface of H 2 /Pd(100) has been determined using the density- 
functional theory together with the generalized gradient approximation (GGA) ]3T| and the 



full-potential linear augmented plane wave method [p2| , |33l . Ab initio total energies have 



been evaluated for more than 250 configurations and have been parametrized in a suitable 



form for the dynamical calculations |[23|| . The substrate atoms are assumed to be fixed since 
due to the large mass mismatch between adsorbate and substrate for H 2 /Pd there is only 
little energy transfer to the substrate phonons. The quantum dynamics is determined in 



a coupled-channel scheme within the concept of the local reflection matrix (LORE) ||3"4|j3"5 
and the inverse local transmission matrix (INTRA) [Q. This very stable method, which 
has been employed before in a high-dimensional study of the adsorption of H 2 /Cu(lll) 
J37J, is closely related to the logarithmic derivative of the solution matrix and thus avoids 
exponentially increasing evanescent waves which cause numerical instabilities. By utilizing 
all symmetries of the hydrogen wave function it has been possible to effectively include up 
to 21,000 channels per total energy in the dynamical calculations. 

III. RESULTS 



Figure [I] presents the results for the sticking probability as a function of the kinetic energy 
of the H 2 beam incident on a Pd(100) surface [p3]| . The dashed curve, which corresponds 



to H 2 molecules initially in the rotational ground state ji = 0, exhibits a strong oscillatory 
structure for low energies. These oscillations are a consequence of the quantum nature of 
the hydrogen beam flT8| , |38|| . They are smoothed out if the initial rotational population and 
the energy spread typical for molecular beam experiments |24j] are taken into account (solid 
line in fig. M. This curve should be compared with the experimental results of Rendulic et 



al. pj]. The theoretical curve agrees quite well with the experimental data. Although no 
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FIG. 1. Sticking probability versus kinetic energy for a H2 beam under normal incidence on a 
Pd(100) surface. Experiment: circles (fromref. [^4|); theory: H2 molecules initially in the rotational 
ground state (dashed line) and with an initial rotational and energy distribution adequate for 
molecular beam experiments (solid line) (from ref. p3| ). 



precursor state exists in the PES and the energy transfer to substrate phonons is not taken 
into account, the initial decrease of the sticking probability with increasing kinetic energy 
is well reproduced. 

The initial decrease results from a dynamical steering effect which had been proposed 
earlier (see, e.g., P5|J25[] ), but not confirmed theoretically. Molecules approaching the surface 
from the gas phase will be attracted to non-activated paths towards dissociative adsorption 
by the potential gradient. The slower the molecules are, the more likely it is that they 
actually follow these attractive paths. By increasing the kinetic energy the time that the 
gradient acts upon the molecules is shortened. More molecules will then hit the repulsive 
part of the potential without being steered to non-activated paths and will be scattered 



back into the gas phase |p3|| . This causes the decrease in the sticking probability. By further 
increasing the kinetic energy the molecules will eventually have enough energy to directly 
cross the barrier which leads to the increase of the sticking probability at higher energies 
(see fig. [I]). In the quantum dynamical coupled-channel description the steering effect is 
reflected by the fact that at low energies more channels are needed in order to get converged 
results than at high energies. This indicates that there is a strong rearrangement between 
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FIG. 2. Orientationally averaged sticking probability versus kinetic energy for different initial 
rotational quantum numbers ji of the incoming molecular beam. The molecular beams are assumed 



to have an energy spread of AE/Ei = 2Av/vi = 0.2 [24] (Ei and Vi are the initial kinetic energy 
and velocity, respectively). 

the different channels at low energies due to the steering. 

The steering occurs in all dynamical degrees of freedom. Therefore by increasing the 
energy of, e.g., the rotational degree of freedom of the hydrogen molecule the steering and 
thus the sticking probability should be diminished. This can already be seen in fig. |l], where 
the sticking probability of the rotationally populated beam is on the average slightly lower 
as compared to molecules in the rotational ground state. This effect is shown in more detail 
in fig. which displays the orientationally averaged sticking probability 
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(1) 



m;=- Ji 



versus initial kinetic energy for jj = 0, . . . , 4. Figure clearly demonstrates that rotational 
motion hinders sticking, especially at low kinetic energies, i.e., the regime where the steering 
effect is operative. The faster the molecules rotate, the more the dissociative adsorption is 
suppressed, because molecules with a high angular momentum will rotate out of a favor- 
able orientation towards adsorption during the dissociation event. This hindering effect of 
rotations becomes smaller, however, at kinetic energies larger than ~0.2 eV, where direct 
activated adsorption is dominant. 



The suppression of the sticking probability by additional rotational motion can actually 
be used to discriminate between the precursor and the steering mechanism. The precursor 
state is usually assumed to be a physisorption state. There are only little directional forces 
for molecules adsorbed in a physisorption state, they can almost freely rotate ||39|| . The 
trapping probability into the physisorption state and thus the sticking probability in the 
precursor model should be almost independent of the initial rotational state, in contrast 
to the steering mechanism. Unfortunately it is not easy to prepare a molecular beam in 
a single quantum state. However, by seeding techniques the translational energy of a H2 
beam can be lowered in a nozzle experiment without changing the rotational population of 
the beam (the translational energy can not be increased since there is no lighter seeding gas 
than H 2 ). In fig. |3| we have plotted the orient at ionally averaged sticking probability versus 
the rotational temperature for different kinetic energies. Experimentally the rotational tem- 
perature of a H_2 beam can not be lower than the corresponding translational temperature 
(a kinetic energy of 200 meV, e.g., corresponds to a nozzle temperature of 1200 K), however, 
theoretically all combinations of kinetic energy and rotational temperature are feasible. For 
kinetic energies below ~40 meV there is a strong dependence of the sticking probability on 
the rotational temperature. By increasing the rotational temperature the sticking probabil- 
ity can be decreased by more than a factor of two at these kinetic energies which should be 
observable in experiment. At large kinetic energies the suppression is less pronounced which 
could already be inferred from fig. 0. 

Interestingly enough, rotational motion seems to suppress sticking in general in the sys- 
tem n 2 /Pd(100). We have checked this for kinetic energies E,i < 0.45 eV and rotational 
quantum numbers ji < 8. Also the observed rotational cooling in desorption of H2/Pd(100) 
p3| , |40[] supports these findings. This situation is different in the system H 2 /Cu(lll) where 



a non-monotonous dependence of the sticking probability on rotational quantum number ji 
has been observed fl,ffl]]: Rotational motion is found to hinder adsorption for low rotational 
states (ji < 4) and enhance adsorption for high rotational states (ji > 4) [Q]. The enhance- 
ment for high j states is related to the elongation of the molecular bond at the barrier 
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FIG. 3. Orientationally averaged sticking probability versus rotational temperature of the in- 
coming beam for different kinetic energies. 

position in the late barrier system H2/Cu(lll) which leads to a decrease of the rotational 



constant and thus to an effectively lowered barrier for high j states |l^ , p^ , p^ , pr5| , p!6| , ^2| . In 



the system H 2 /Pd(100) these late barriers, however, are absent [21]. 

There is still an effect that can over-compensate for the suppression of the sticking proba- 
bility by rotational motion, namely the orientational or steric effect [23|| . The most favorable 
orientation to adsorption is with the molecular axis parallel to the surface. Molecules with 
azimuthal quantum number m = j have their axis preferentially oriented parallel to the 
surface. These molecules rotating in the so-called helicopter fashion dissociate more easily 
than molecules rotating in the cartwheel fashion (m = 0) with their rotational axis preferen- 
tially parallel to the surface since the latter have a high probability hitting the surface in an 
upright orientation in which they cannot dissociate. This steric effect, which has also been 



investigated in a number of model studies for purely activated adsorption ||ll|ji^jl^ji^Jlq . [16 



can clearly be seen in fig. |] where the sticking probability for one fixed kinetic energy of 
Ei = 0.175 meV is plotted. Indeed the rrii = ji data even rise with increasing quantum 
number ji at this relatively high kinetic energy, while the m, = and the orientationally 
averaged results are decreasing. At lower kinetic energies (which is not explicitly shown 
here) , where the steering is more pronounced, also the rrii = ji data decrease. 
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FIG. 4. Sticking probability versus initial rotational quantum state j{. Diamonds: orientation- 
ally averaged sticking probability (eq. [I]), triangles: rrtj = (cartwheel rotation), circles: rrtj = j% 
(helicopter rotation). The initial kinetic energy is E{ = 0.175 eV. 

IV. CONCLUSIONS 

In conclusion, we have reported a six- dimensional quantum dynamical study of adsorp- 
tion and desorption in the system H2/Pd(100) using an ab initio potential energy surface. 
We have shown that the initial decrease of the sticking probability with increasing kinetic 
energy is due to dynamical steering. We have focused on the steering effect in the rotational 
degree of freedom of the hydrogen molecule and shown how the steering effect can be further 
confirmed experimentally. Our study demonstrates that the combination of ab initio poten- 
tial energy surfaces with high-dimensional quantum dynamical calculations can lead, due to 
the microscopic information, to a quantitative as well as new qualitative understanding of 
processes at surfaces. 
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